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SUMMARY 

In this report a formula for calculating the induced drag of multiplanes with end plates is 
derived. The frictional drag of the end plates is also calculated approximately. It is shown fft&t 
the reduction of the induced drag, when end plates are wsedj is sufficiently large to increase the 
efficiency of the wing. 

Curves shovnng the reduction of drag for monoplanes and biplanes are constructed; the influence 
of gap-chord ratio } aspect ratio, and height of end plate are determined for typical cases. The 
method of obtaining the reduction of drag for a multiplane is described. 

Comparisons are made of calculated and experimental results obtained in wind tunnel tests 
with airfoils of various aspect ratios and end plates of various sizes. The agreement between cal- 
culated and experimental results is good. 

Analysis of the experimental results shows that the shape and section of the end plates are 

important. ^ 

INTRODUCTION 

The end plates which are dealt with in this report are fins or shields which are attached 
to the tips of airfoils; the plane of the fin is perpendicular to the span of the airfoil. The purpose 
of this device is to serve as a barrier to the flow along the span and around the tips of the airfoil. 
This flow in a vertical plane containing the span is usually called the transverse flow. By 
obstructing this transverse flow its kinetic energy is diminished and as a consequence the 
induced drag of the airfoil is thereby reduced. 

From wind tunnel tests made at Gottingen (Reference 1) and at the Langley Memorial 
Aeronautical Laboratory (Reference 2) ? it has been learned that the induced drag may be 
materially reduced by the use of end plates. The reduction of the induced drag exceeds the 
additional frictional drag of the end plates which were used in these tests for all but small 
values of the lift coefficients. 

The German tests were analyzed but the method of analysis gives no direct information 
on the value of the frictional drag of the end plates. The American tests have not been pre- 
viously analyzed. 

In this report a formula for calculating the induced drag of multiplanes with end plates is 
derived. The induced drag is calculated by finding the kinetic en'ergy of the transverse flow 
using the method of conformal transformation. The induced drags of the monoplane and 
biplane with end plates are found as special cases of this general formula. 

The formula for the monoplane is used to calculate the reductions of drag for the airfoils 
tested at ♦Gottingen and the Langley Memorial Aeronautical Laboratory. Comparison of 
calculated and experimental results for airfoils of various aspect ratios with end plates of various 
sizes are made. Curves illustrating the total reduction of drag for monoplanes and biplanes 
are constructed; the influence of such major factors as gap-chord ratio, aspect ratio, and height 
of end plates is determined. In all cases the frictional drag of the end plates is considered in 
calculating the total drag. 
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INDUCED DRAG OF A MULTIPLANE WITH END PLATES 

The longitudinal projection of the multiplane consists of equal, parallel straight- lines. 
The same projection of the end plates will also be considered as straight lines whose directions 
are at right^ angles to the projection of the multiplane elements. 

The induced angle of attack and consequently the induced drag depend upon the induced 
downwash. It is sufficient- then for the purpose in hand to consider the transverse flow in a 
vertical plane at right angles to the air stream. 

This two-dimensional transverse flow may be determined in several ways. In this report 
it has been obtained by using the method of conformal transformation. The rectilinear polygon 
bounding the longitudinal projection of the wings and end plates is transformed into a single 
straight line. The flow around such a line is well known. 

Let the plane of the lines with end pieces be thejplane of a complex variable z as shown 
in Figure 1. The mid-point of the system of lines is at the origin and the direction of the lines 

coincides with that of the axis of imagi- 
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naries. The end pieces CE and C'E' 
are then parallel to the real axis and are 
each of length 2Ji. The bounding rec- 
tangle BFF'B' has dimensions 2b'X2h'. 
It is assumed that h^Ji'. 
" The transformation to~an auxiliary 
f-plan-e is effected by means of the 
theorem usually called the Schwarz- 
Christoffel theorem. According to this 
theorem, if z=x+iy and f = £ + t7? then 
any polygon bounded by straight lines 
in the 2-plane can be transformed into 
the axis of £, points inside the polygon corresponding to points on one side of the axis of f. 
The transformation is 
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Fig. 1 — Conformal representation of biplane with end plates 
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(1) 



where a h a g , * • • - a n are internal angles of the polygon in the 2-plane and % u f f? - . . £ tt are 
the points on the axis of £ that -correspond to the angular points of the polygon in the z-plane 
and A is a real constant. 

In our application we shall consider the polygon -coABCDEF-f od in the 2-plane as 
transformed into the upper half of the f-plane. By the theory if the point ?=od be taken to 
correspond with one corner of the polygon the corresponding factor in ..the expression for dz/dt 
is omitted. Also we may choose three values of £ arbitrarily since a polygon similar to a given 
polygon of n sides may be constructed by means of n— 3 relations between the sides. 

We shall accordingly choose the three points as follows: 



D F 
o 1 



Then to the points E and G will correspond the values f=c and £=p respectively. The 
transformation is then found in the form 



dz 

di = a 



where a is a real constant. 



(2) 
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We may integrate (2) by substituting t=mu (mod. k). snco denotes an elliptic function 
and from the relations given in textbooks on these functions we find, 



,=|[z( w )-c.(l-^i^)]+iS 



(3) 



Here Z (w) is Jacobi's zeta-function. 

The correspondence between the various points is exhibited in the following table: 

A B CDE FGco 

z-¥ 7 -V+il't -h + ib', h + ib', W+iV, h' 7 

t,-i/h -1, -c, 0, ' c, 1, 

a> } -E+iE', -E, -a 0 0, a 93 E, E+iE' s iE' 

E and E are the complete elliptic integrals of the first and second kind, E f being the complete 
elliptic integral of the first kind for the comodulus ¥ (i' 2 = 1— i*}. 

E 

Since c = sn u 0 we may now use the relation 1— y — V sn s w 0 = Z' (ca 0 ). This relation is 

* 

also found in textbooks on elliptic functions. The prime accent indicates differentiation with 
respect to the variable We then have 

z== Te [~°° Z ' M 0)1 + & * and j8 constants. (4) 
Using (4) and the correspondence given in the table we have the following equations: 

j [~a> 0 Z' O,) +Z (%)] + |8 = a+i5' (5) 
£l-KZ' M] + (3=V+ib' (6) 
| [- {K+iK')Z> -i* + j8«i' - (7) 

[Z (e^-^Z' (<*)]; i&'-fl 
V=-fl(E) Z' M] 



Solving these we have: 



h_ _ Z (<j]q) — Cd a Z' ( cj 0 ) 

¥ K'Z'{u 0 ) + - c 
7i'_ -KZ' (gg) 



If B /-Y, Z' ^) = 2 g(Z+ T Z 7 ) (10) 

(11) 

.-V ^V^ ZM+lg+i] ' (12) 



a 2¥ (K+ K'y) 
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The transverse flow around the lines with end pieces is parallel to the real axis at infinity 
in the 2-pIane, the direction of flow at infinity being in. the direction of the positive cc-axis. The 
potential function for this flow is more easily expressedin terms of t. It is $ 1 = <p t + ifa = a Vt } 
where V is the velocity of the flow in the s-plane at a great distance from the origin. 

If we superimpose on this flow one whose potential function is — Tz we shall have the 
potential function for the flow set up by the motion of the lines in'a fluid otherwise at rest. 
The lines will move in the direction of the positive real axis in the 2-plane with a constant 
velocity V. 

To find the kinetic energy we evaluate the integral 

— 2 pj P ^~ ds (Reference 3, p. 44). 

The path of integration, L } is the complete boundary of the wing, jp is the real part of (a Vt— Vz) i 
gjj is the normal velocity component along the boundary and ds is the element of length of the 
boundary. We have 

t = 2 '™![K 1+ p)- c *] 

TP 

By using again the relation 1 -y— ¥c* = Z' as well as equations (10) and (11), we have 

1 2 p[ tt IK L E(K+yK'U { } 

This may be written T= ^ p V*A M where 

A M - L I ''KXE+WTU (U) 

Prandtl (Reference 4, p. 48) has shown that^if it is assumed that all air particles within an air 
stream of cross sectional area A M are uniformly deflected downward by the wing and all outside 
are not- deflected, then the correct" lift and the correct work due to the induced drag are obtained 
by application of the impulse theorem and the energy theorem. Doctor Zahm has suggested 
that A M be called the area of the equivalent air stream and this term will be used throughout 
the report. 

As the minimum induced drag of the wing system is known to be D^- A — -t— j where L is 

4g Am 

the lift and q the dynamic pressure, it may be computed by using equation (14). 

THE MONOPLANE WITH END PLATES 

The area of the equivalent air stream of the monoplane with end plates is found from 
(14) by setting y==0. 

K*h*r2E 1 

If the span is b we have A M =—^-\ (15) 

This reduces when Ji r =h = 0 to 

A M =^-~ which is known to be correct. 

The ratio of the area of the equivalent air stream of the monoplane to that of the mono- 
plane with end plates is 

# = /op \ ' - -(16) 



J= g+ 0.6515 (jf) is also given in Table I. 
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Values of R were calculated in the following manner: When 7 = 0, Z f (co 0 ) =0, and a 0 was then 
found from a table of elliptic functions (Reference 5) after the modulus Tc had been arbitrarily 
assigned. Hs usually expressed as the sine of a modular angle, 0. Then 2h[b and R were calcu- 
lated using equations (8) and (16). 

Table I contains the values of R for various values of 2h/b. The quantity / as given by 
Nagel (Reference 2) is related to R by the equation 

The value of R obtained by using the approximate formula 

Let us consider now an airfoil whose span is h and chord c. We shall assume that the longi- 
tudinal projection of the airfoil is a straight line. The width of the end plate is taken as a chord 
length where no statement to the contrary is made. The coefficient of minimum induced drag 
is then 

(7i=lift coefficient 
S = area of wing. 

Using (16) this may be written 

The use of end plates introduces an additional frictional drag. This frictional drag may be 
designated D F and is given by 

B F =0 F qS r 

The end plates act as airfoils so that G? is a drag coefficient and S' is the area of the two end 
plates. 

The difference of total drag coefficient, C Df and profile drag coefficient, C Dpt is 

n _p _0 L *RS n S' 
In O^-— . b r -g 

C £ S 

The same expression for a wing without end plates is To determine the actual gain with 

a given set of end plates we subtract these two expressions. If we call this difference AC D , we 
have 

AC fl = ^f (!-£)- (18) 

If the width of the end plate is a chord length and the height 2ft, then 

^_4cft_4ft 
S ~ cb~ b 



Then AC D =^(l-5)-2^(f) 

rv, 



(19) 



If we use i?=fTT= 



we have 



^L-*C,(f) (20) 



(¥) 
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From this the reduction of drag for a given end plate may be easily found. Figure 2 shows 
parabolas obtained from (20) by using a fixed value of b~/S } C F} and different sizes of end plates. 
For smaller values of b 2 fS the slopes of the parabolas increase but the intersection points on 
the vertical axis are not altered. There is consequently a greater average reduction of drag for 
the smaller values of h 2 /S o.ver the usual range of values of the lift coefficients 

ANALYSIS OF WIND TUNNEL RESULTS 

The Gottingen tests (Reference 1) were made with single airfoils. The end plates were 
flat and of various shapes and sizes. Two aspect ratios were used, 8/3 and 4/3. In the analysis 
made of these tests the frictional drag of the end plates was not actually calculated. In this 
report this quantity has been calculated approximately. 

Figures 3, 4> and 5 show the results of the Gottingen tests, for aspect ratio 8/3, in the 
form of polar curves. In each case except for end plate 3, Figure 4 7 the reduction of drag was 
computed from formula (19) using the appropriate value of aspect ratio, ratio of length of end 
plate to span and area of end plates. In the case of end plate 3, since the end of the wing was 
only partially shielded, an average height of the plate, over the chord, was used. The value 
of C F which gave the best agreement for the frictional drag of the end plates is 0.0140. 

This calculated reduction of drag was then 
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laid off on each diagram starting from the polar 
obtained when using the whig without end 
plates. Unfortunately, no tabular values of 
drag and lift coefficients were given in the 
Gottingen report. Consequently, no accurate 
comparisons of calculated and observed reduc- 
tion of drag could be made. 

On the who!e ; the agreement between cal- 
culated and observed reductions of drag is 
good for that portion of the polar curve where 
the angles of attack are smaller than the 
burbling angle. The use of a partial end plate 
such as end plate 3 is not satisfactory. In this 
instance the disturbance of flow around the end 
of the wing apparently increases the frictional 
drag at low values of C L and the reduction of 
drag is not" as great at the higher values of C L . 

The American tests (Reference 2) were 
made at the Langley Memorial Aeronautical 
Single airfoils, wing section N. A. C. A. 73, with aspect ratio 6 were used in 
Only two plates, one circular, the other trapezoidal, were tried. These plates were 
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Fig. 2.— Monoplane with end plates. Calculated reduction of drag, 
ACz>, against lift coefficient for various values of 2kfb. 2?t=neight of 
end plate. &=span. Aspect ratio, 6. Drag coefficient, O=G.01 



Laboratory, 
these tests.' 

streamlined and consequently the value of C F to be used in making the calculations is different. 
The value 0.008 (Reference 1) proved satisfactory. 

The results of these tests are shown graphically in Figures 6 and 7. The reduction of drag 
was calculated as before, C F only being changed in value. Since tabular values are available 
much closer comparison could be made. The calculated and observed reductions of drag are 
shown in tabular form in Table II, and graphically in Figures 6, 7, and 8. The agreement is 
good, considering that we have only an approximation to the acutal reduction in consequence 
of the simplifying assumptions made in the calculations. 

Comparing the two series of tests we see that theory and experiment show a larger reduc- 
tion of drag for smaller aspect ratios. The use of end plates with good aerodynamical profile 
sections is preferable since this reduces the drag for all values of 0 L , This is illustrated by the 
two series of tests studied; in one the coefficient is 0.014, and in the other 0.008. 

The agreement between calculated and observed reductions of drag in the American tests 
suggests using an end plate of given area in the best possible way. Since the use of the entire 
height of the trapezoidal end plate (fig. 6) in the formula gave calculated results agreeing very 
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Fig. 3 —Polar curves of Gottingen tests showing reduction of drag 
due to end plates 2 and 3. Aspect ratio 8/3 



Fig. 6.— Polar curves of N. A. C. A. tests showing reductions of drag due to disk end 
plates. N. A. C. A. 73 airfoil. Aspect ratio, 6 




Fig. 4 —Polar curves of Gottingen tests showing reduction of drag 
due to end plates 3 and 7. Aspect ratio, S/3 
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Fig. 7— Polar curves of N. A. C. A. tests showing reductions of drag due to trapezoidal 
end plates. N. A. C. A. 73 airfoil. Aspect ratio, 6 
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Fig. 5.— Polar curves of Gottingen tests showing reduction of Fig. S.— Calculated and observed differences of drag of wings with different end plates, 
drag due to end plates 4 and 6. Aspect ratio, S/3 N. A. C, A. 73 airfoil. Aspect ratio, 6 
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Fig. II. — Biplanes with end plates. Calculated reduction of drag, LCd, plotted 
Fig. 9. — Calculated redaction of dragging disk end plates and increasing wing area against lift coefficient for various values of k/G. 2ft=lieigbt of end plate. 2<7— 
by an amount eqaal to area of end plates. N. A. C. A. 73 airfoil. Aspect ratio, 6 gap. Span/chord = 6. Gap/chord=I.0Q. Drag coefficient, CV*=0.01 



AG, 



.06 
.05 
.04 
.03 
.02 
.Of 
0 

-or 























































h/i 








4 




































/ 


2 
/ 




































































f 





































































0.2 0.4 0.6 



0.8 f.O 



.05 



-.04 



.03 



.02 

.Of 
0 

-Of 



1.2 ... f.4 • L6 



-.02- 
— 0 























/ 
























} 




I 
































2/ 


























































/ 

























































































0.2 _ 0.4 as 



{.2 _ f.4 



Fig. 10— Biplanes with end plates. Calculated reduction of drag, & Cd, plotted 
against lift coefficient for various values of UfG. 2ft=h6igrit of end plate. 2<?= 
gap. Span/chord =6. Gap/chord=0.6. Drag, coefficient, 0=0.01 
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Fig. 12.— Biplanes with end plates. Calculated reduction of drag, ACd, 
plotted against lift coefficient for various values of hfO. 2A- height of 
end plate. 2£=gap. Span/chord =6. Gap/chord~L40. Drag coeffi- 
cient, 0=0.01 



well with observed results, it appears that an end plate of a given area, and therefore having 
a fixed Motional drag, may have a greater effect on the drag in one shape than in another. 
In particular, the tests and calculations with this end plate make it seem plausible that the 
greater part of the restriction of the transverse flow occurs over the forward portion of the 
span and that a shape which has greater height in the portion nearer the leading edge is advan- 
tageous to-use. - Zl. . 

Nagel (Reference. 1) has stated that the reduction of drag when using end plates, does not 
equal the corresponding reduction obtained when thgLspan is increased by an amount equal to 
the sum of the heights of the end plates. The reason for adding to the span an amount equal 
to the total height of end plates is not quite evident. The addition to the span might also be 
made, with as good reason, by adding to the wing area the sum of the areas of the two end 
plates. Such a comparison is. shown graphically in Figure 9. The polar curve of an airfoil is 
shown compared with the computed polar curves when end plates aA added to the wing and 
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when the wing area is increased by an amount equal to the total area of the end plates. The 
comparison was made using the same drag coefficient for the end plates as for the additional 
wing area. The results do not differ much. 

CALCULATED REDUCTION OF DRAG FOR MULTIPLANES WITH END PLATES 

The method will be used to calculate first the reduction of drag for a biplane. The area of 
the equivalent air stream for a biplane with end plates is given by equation (14). The end 
plates must have heights equal to or greater than the gap of the biplane, for otherwise the flow 
would change to such an extent that (14) would no longer hold. This formula reduces when 
A = fi/, that is, when the end plate has a height equal to the gap of the biplane, to one which has 
been previously derived (Reference 6). It is evident from the correspondence of values between 
t and co that when 

l~V, a 0 =K } Z(a 0 )=0 and Z' (6>*) = 1-J-F- 

From (10) we find then 

E-El rs A A „ ¥ s 
^ E'-K'V and = wb * E'-VK'* 
For a span 6 this reduces to 

a 

Am ~ 4 ' E'-V E' 

which is the value previously found. 

The expression for reduction of drag is derived in a manner quite similar to that used for 
the monoplane. We find 



& = span of biplane. 
c = chord. 
0 L = lift coefficient. 
.4^ = area of the equivalent air stream of a biplane without end plates. 
Ax =area of the equivalent air stream of a biplane with end plates. 
CV = drag coefficient of end plates. 
2h = height of end plate (width assumed to be c). 

If we introduce PrandtPs ratio x ^^jT and a s * m ^ ar ra ^° 



Xi = - A — i — we have 

A ^ = ^|(x-^-^(^) (22) 

The value of x x depends on both the gap/chord ratio and the ration of height of end plate to gap. 
The former ratio is y and has been chosen in each calculation. The calculations were made in 
the following manner: 

XT 
9 

(a) 7/ (osq) was found from 7J («„) = K (E^y E') 7 after 7 tas beeU arbitraril F assigned. 

E 

(b) From dft'c^-^-hZ' (co 0 ), u 0 was found using the table of elliptic functions (Reference 4). 

(c) The value of f = Z ^~^ 7/ (cJg) was calculated next. 

b ttk'+E'Z' (co,) 

(d) A Ml - — K ' y ^ [ 2 jf ~ ~ 2 7/ fa)] was evaluated then quite easily. 

Figures 10, 11, and 12 show the calculated reductions of drag for span/chord of 6 and 
gap/chords of 0.6, 1.0; and 1.4. In each case the following values of Ti[hf were used; lt(k f = 
1, 2, 4, 6 and, as in Figure 2, the value G F — .01 was used as it represents an average value of C F . 
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The curves illustrate how properties of end plates affect biplanes of various proportions, 
The increase of drag at lower values, of C L is considerably less than it is for larger values o£ 
gap/chord while the reduction of drag, especially when larger plates are used, is relatively large. 

The reduction oj drag persists at the higher values of Q L for plates which are quite large, 
so large as to be useless in practice. 

The method used in calculating the effect of end plates for biplanes may be used for multi- 
planes. Formula (22) would differ only as far as x is. concerned. This quantity may be cal- 
culated, however (Reference 7), and the reductions of drag for multiplanes may be calculated 
in precisely the same manner as that just described and carried out for biplanes. 

CONCLUSIONS 

Calculations show that the induced drag of monoplanes and multiplanes may be decreased 
by attaching end plates to the ends of the wing. The frictional drag of the end plates may be 
calculated approximately. The reduction of the induced drag exceeds the additional frictional 
drag due to the end. plates at all but small values of the lift. For given dimensions of wings 
and end plates the reduction of drag less the frictional drag of the end plates varies directly 
as the square of the absolute lift coefficient. . The average reduction of drag increases, as the 
aspect ratio decreases. Calculations and experiments agree quite satisfactorily for single wings 
equipped with end plates. 

Wind tunnel tests show that the coefficient used, in calculating the frictional drag of the 
end plates may be reduced materially by fairing the end plates. The shape of the end plate 
determines to some, extent the reduction of induced drag. Further tests would be necessary 
to find end plates of good shape and having at the same time a low drag coefficient. 

Recent experiments have shown that much higher lift coefficients can be obtained than 
has been the case up to now with the conventional airfoils. Since the reduction of drag due to 
end plates is much greater at the higher values of the lift coefficient^ the use of end plates 
together with the use of means for increasing the lift will result in a material improvement of 
the efficiency of the present-day airplane. 

Langley Memorial Aeronautical Laboratory, 

National ^Advisory Committee for Aeronautics, 
Langley Field, Va., January 20, 1927. 
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TABLE I 

Values of R and 2hfb for given values of 6. 

#=(area of the equivalent air stream of a monoplane) /(area of the equivalent air stream of a monoplane 
with end plates) . 
2h/b= (height of end plate) /(span of wing). 



$ degrees 


2hlb 


B 


JNaget s 
formula 


.Per cent 
difference 


15 


0. 0173 


0. 971 


0.974 


+G. 33 


on 


. Uoli 


. yio 


. yoi 




25 


.0493 


. 909 


.925 


L 76 


30 


. 0721 


. 877 


.895 


2. 05 


35 


.0996 


. 844 


. 859 


L 78 


40 


. 133 


.soo 


. 820 


2. 50 


45 


. 173 


.751 


.777 


3. 47 


50 


.222 


.704 


.731 


3. 83 


55 


.285 


. 653 


.678 


3. 80 


60 


.349 


. 610 


.633 


3. 7 7 


65 


.433 


.562 


.582 


3.56 


70 


. 541 


.508 


.527 


3. 74 


j 75 


. 686 


. 452 


.468 


3. 15 


80 - 


. 897 


. 395 


.401 


-j- L 52 


85 


1. 27 


.323 


.321 


-.62 


! 88 


L 78 


.261 


. 252 


-3. 45 


89 


2. 17 


.227 


.217 


-4 40 



TABLE II 

OBSERVED AND CALCULATED REDUCTION OF DRAG 



[Tests made at Langley Memorial Aeronautical Laboratory. Monoplane, 4 by 24 inches; wing 

section, N. A. C. A. 731 



C L 


Cb (disk end plates) j 


Difference of 
observed and 

calculated 
reduction as 

per cent of 
drag 


Cd (trapezoidal end plates) 


Difference of 
observed and 

calculated 
reduction as 

per cent of 
drag 


Observed 


Calculated 


Observed 


Calculated 


0 


' -0.0026 


-0. 00210 


-1-3. 90 


-0. 0010 


-0.00113 


-0. 18 


. 1 


-. 0020 


— . 00196 


+.35 


— . 0007 


— . 00105 


— 3. 50 


. 2 


-.0013 


-. 00159 


t— 2. 40 


— . 0007 


— . ooos 


-.83 


. 3 


0003 


— . 00095 


-4. 60 


— . 0003 


— . 00039 


— . 63 


< 4 


+. 0005 


-. 00005 


+ 3. 15 


+. 00020 


— . 00018 


+ - 11 


! 5 


+.0012 


+. 00110 


+. 95 


.0007 


+ . 00092 


— L 00 


. 6 


• +. 0023 


-K 00250 


- 72 


+. 0016 


+. 00183 


-. 82 


. 7 


-K 0038 


. 00412 


— . 90 


+. 0027 


+. 00289 


— . 54 


. 8 


-K 0064 


.00604 


+'.82 


+. 0036 


+. 00413 


— L21 


.9 


+. 0093 


. 00826 


+ 1. 82 


+. 0060 


+. 00553 


+. 83 


1.0 


+. 0180 


. 01060 


+ S.22 


+.0162 


+. 0070S 


+ 10. 10 



